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This paper describes the pre- flight stage separation analysis that was conducted in sup- 
port of the final flight of the X-43A. In that flight, which occurred less than eight months af- 
ter the successful Mach 7 flight, the X-43A Research Vehicle attained a peak speed of Mach 
9.6. Details are provided on how the lessons learned from the Mach 7 flight affected separ a- 
tion modeling and how adjustments were made to account for the increased flight Mach 
number. Also, the procedure for defining the feedback loop closure and feed-forward pa- 
rameters employed in the separation control logic are described, and their effect on separ a- 
tion performance is explained. In addition, the range and nominal values of these parame- 
ters, which were included in the Mission Data Load, are presented. Once updates were 
made, the nominal pre- flight trajectory and Monte Carlo statistical results were determined 
and stress tests were performed to ensure system robustness. During flight the vehicle per- 
formed within the uncertainty bounds predicted in the pre -flight analysis and ultimately set 
the world record for airbreathing powered flight. 




Nomenclature 

a 

= 

Angle of attack, deg 

P 

= 

Sideslip angle, deg 

ftmax 

= 

Peak sideslip angle, deg 

C A 

= 

Axial force coefficient 

C N 

= 

Normal force coefficient 

C m 

= 

Pitching moment coefficient 

§elev 

= 

Elevator deflection, deg 

&rud 

= 

Rudder deflection, deg 

da/dt 

= 

Rate of change of angle -of-attack, deg/s 

♦ 

= 

Body roll angle, deg 

Fx 

= 

Axial Force in vehicle body frame, lbf 

P 

= 

Body roll rate, deg/s 

q 

= 

Body pitch rate, deg/s 

r 

= 

Body yaw rate, deg/s 

a 

= 

standard deviation 

XSEP 

= 

Horizontal separation distance, inches 
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I. Introduction 


O n November 16, 2004 NASA’s X-43A demonstrator vehicle, the first free flying air-breathing scramjet pow- 
ered aircraft, successfully completed its final test flight, reaching Mach 9.6 - a world speed record for powered 
flight. This flight occurred less than 8 months after the X-43A’s first successful flight, in which the planned peak 
speed of Mach 7 was attained. 1 " 2 

One of the highest risk events during both of these flights was the separation of the X-43 A research vehicle (RV) 
from the Hyper- X Launch Vehicle (HXLV). Although stage separation of two non-axisymmetric bodies at these 
flight conditions (Mach 7 and 10, dynamic pressure of 1000 psf) had never been done, in both instances stage sepa- 
ration was successful and the RV reached its target attitude prior to the start of the test sequence without re- 
contacting the HXLV. 



Figure 1. Artist concept of X-43 A stage separation. 


Even though the success of the stage separation in the first X-43A Mach 7 flight increased the overall confidence 
in the stage separation simulation modeling, there were some minor differences that were observed between the 
flight data and pre -flight models. Efforts were made to understand these differences by systematically reconciling 
the simulation models with the flight data and incorporating the relevant findings and model adjustments into the 
Mach 10 separation simulation. In addition, these flight data comparisons were used to reassess model uncertainty 
levels and define additional stress cases that were run to ensure system robustness. Finally, since the vehicle had less 
stability and control effectiveness at Mach 10, both Computational Fluid Dynamics (CFD) and wind tunnel results 
were used to modify the aerodynamic interference model to better match the Mach 10 flight conditions. The changes 
in the vehicle aerodynamics were significant enough to require that the separation control logic be redesigned. This 
paper describes the development of the overall Mach 10 separation strategy and the corresponding risk reduction 
activities that were undertaken to account for the different flight condition and demonstrate that the system was ro- 
bust to any errors seen in the previous Mach 7 flight. 

II. Mach 10 Separation Simulation Development 

A. Simulation Methodology 

To assess the risk associated with stage separation, a 14 degree-of- freedom simulation had been created early in 
the X-43 A program using the Automatic Dynamic Analysis of Mechanics Systems (ADAMS) Simulation software. 3 
This simulation was composed of various user-defined modules that modeled the vehicle aerodynamics, ejector pis- 
tons, control system, inertial navigation system, and control surface actuators. The entire simulation, which includes 
the user-defined modules wrapped around the ADAMS solver, is known as SepSim. 7 In addition, another independ- 
ent separation simulation was developed prior to the Mach 7 flight as an additional risk reduction measure to verify 
the SepSim results. This independent simulation, which was developed using the Program to Optimize Simulated 
Trajectories (POST) 4,5 , was less detailed than SepSim, but was simpler, more concise and ran faster, while providing 
similar results. 6 

In order to take advantage of the strengths of each of these simulations, the Mach 10 updates were applied to 
both POST and SepSim. Because of its rapid Monte Carlo capability and efficiency as a parametric analysis tool, 
POST was used as the “work-horse” simulation while the higher fidelity SepSim, which, unlike the POST simula- 
tion, included high-order actuator and sensor models, served as the final truth model. Each simulation was verified 


2 

American Institute of Aeronautics and Astronautics 




separately using standalone models. In addition, both nominal and off-nominal cases were run in both POST and 
SepSim and the results were compared and shown to be in good agreement. Similar comparisons were made peri- 
odically throughout the Mach 10 study to ensure that traceability between the two simulations was maintained. 

B. Mach 10 Separation Timeline 

The entire Mach 10 stage separation event was 2.5 sec in duration. The nominal separation scenario is shown in 
Fig. 2. The event begins as two pyrotechnically actuated pistons push the vehicles apart. At 100 ms the pistons reach 
the end of their 9 -inch extensions and the HXLV commands a 10 deg elevon deflection to induce a nose-down pitch- 
ing motion, moving it further away from the separated RV as it descends uncontrolled back to Earth. For roughly the 
first 250 ms the two vehicles are in close enough proximity that aerodynamic interference effects have a significant 
effect on the trajectory. These effects gradually diminish and are gone by 500 ms when the horizontal separation 
distance between the vehicles is approximately 110 in. Figure 2 shows the times of the various loop closures with 
rate feedback loops closing at the end of the piston push and outer loops (a and 9) closing at the beginning of free 
flight. The yaw rate loop is faded in to improve off-nominal performance by balancing the likelihood of sideslip 
excursions and high rudder actuator rates. 

The primary objective was to maximize the probability of a successful separation. There were three specific re- 
quirements that were the primary means of defining separation success: 1) to avoid a re -contact between the two 
flight vehicles, 2) to minimize the likelihood of a loss of control resulting in large attitude excursions (i.e., a or (3 
exceeding +10°), and 3) to reach the attitude target box required for the engine test at the end of the 2.5 sec separa- 
tion event (a = 1.0° + 1.0°, [3 = 0.0°+ 0.5°) with zero attitude rates. 



0 msec 

Start of Separation 

Mach = 9.6 

Dynamic pressure = 1000 psf 
a = (3 = 9= 0 deg 

100 msec 

End of Piston Push (end of 9 in. extension) 
Pitch & roll rate feedback loops closed 
Begin fade-in of yaw rate feedback loop 
HXLV elevons to 1 0 deg at 1 20 msec 

250 msec 

Beginning Transition from Interference to 
Free Flight 

Xsep = -44 in. 

500 msec 

Free Flight 

Xsep = -1 10 in. 

Outer feedback control loops closed 

900 msec 

Yaw rate feedback loop closed 

2500 

msec 

End of Separation/Start of test 

a = 1 .0 deg, (3 = 9=0 deg 


Figure 2. Nominal Mach 10 separation ti meline. 


C. Simulation Model Updates 

In general, the structure of the Mach 10 stage separation simulation was the same as that used to conduct the pre- 
flight analysis for the Mach 7 flight; however, a number of updates had to be made to account for the different vehi- 
cles and flight condition. In addition, model reconciliation studies were conducted in which key simulation models 
were compared with Mach 7 flight data to assess how well the flight data was captured by the pre -flight Monte 
Carlo uncertainty bounds and, if necessary, determine how to adjust model parameters to best fit the observed re- 
sponse. 
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The separation conditions shown in Table 1 were similar to the Mach 7 conditions with the exception of the 
separation Mach number. As with the Mach 7 pre -flight Monte Carlo analysis, extremely conservative uncertainty 
models were assumed. For instance, even though the HXLV was predicted to deliver the RV to a separation angle - 
of-attack between -0.6 and 1.0 deg, this uncertainty was assumed to be as high as +1.5 deg. A similar conservative 
approach was taken with the other initial separation condition uncertainties listed in Table 1. 


Table 1. Mach 10 separation conditions and uncertainty. 


Separation Condition 

Mean 

+3 -a 

Dynamic Pressure, psf 

1066 

210 

Mach Number 

9.6 

0.6 

Flight Path Angle, deg 

1.5 

1.0 

Angle -of-attack, deg 

0.0 

1.5 

Sideslip Angle, deg 

0.0 

1.0 

Body Roll Angle, deg 

0.0 

15.0 

Body Pitch Rate, deg/s 

-0.18 

3.0 

Body Yaw Rate, deg/s 

0.0 

3.0 

Body Roll Rate, deg/s 

0.0 

6.0 


A number of mass property updates were required to account for differences in both vehicles with respect to the 
Mach 7 flight. The largest change was due to the fact that the HXLV did not need to carry ballast to reach the Mach 
10 separation conditions. Consequently, at separation the HXLV mass was lower and there was a significant shift in 
eg position. The nominal weight of the entire HXLV/RV stack at the initiation of stage separation was 10,415 lb 
with the RV alone comprising 2,823 lb of the total. 

The locations and force profiles of the ejector pistons were the same as in the Mach 7 flight. As described in 
Reference 6, the separation process was initiated by two pyrotechnically activated ejector pistons that emanated 
from the HXLV and remained in contact with the RV for approximately 0.1 sec until reaching the end of their 9-in. 
stroke length. The pistons were modeled in the simulation as a spring with a test-derived axial force profile that was 
a function of time and extension length, and with lateral restoring forces computed from beam deflection theory. 7 
After the Mach 7 flight the piston model was assessed by comparing the axial acceleration from simulation results to 
acceleration data derived from telemetry. The flight data was well captured in the simulation results as illustrated in 
Fig. 3. Additional comparisons were performed that indicated the piston extension time histories were also consis- 
tent with simulation predictions. As a result of this agreement, no modifications were made to the simulation piston 
models. 

The aerodynamic model employed in the 
separation simulation was a complex database 
constructed from wind tunnel data and comp u- 
tational fluid dynamic calculations that pro- 
vided force and moment coefficients for three 
different vehicle configurations (stack, HXLV 
and RV) in regions of free-flight and aerody- 
namic interference. 8 " 10 Two different eight- 
dimensional databases were used to compute 
the interference aerodynamics on each vehicle 
depending on the value of X S ep (axial separa- 
tion distance). The vehicles were assumed to be 
completely in the interference region from 
separation until Xsep = -44 in (note that Xsep is 
initially 0 when separation commences and 
becomes more negative as the vehicles move 
apart). In the Mach 7 analysis, when X S ep 
reached -69 inches the vehicles were assumed 
to be out of the interference region and a differ- 
ent free flight database was used to compute the 



Time, sec 

Figure 3. Mach 7 axial acceleration profile during piston 
firing, flight data vs. pre-flight Monte Carlo. 
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aerodynamic forces and moments. In the region between -44 and -69 inches, the base coefficients were ramped line- 
arly from the interference values to the free flight values. The aerodynamic coefficient build-up in the simulation 
also included Monte Carlo uncertainty terms to account for potential errors in the model. In the interference region, 
the uncertainty of each coefficient was modeled as an additive increment that was a function of X S ep- In the free 
flight region, a multiplier was applied to each nominal coefficient and an incremental bias term was also added. 

Even though the Mach 7 flight had been conducted without any major problems during separation and the RV 
achieved an acceptable attitude prior to the initiation of the engine test sequence, the vehicle transitioned to the com- 
manded angle -of-attack slower than what had been predicted in pre -flight simulations (see Fig. 4). Moreover, the 
response was outside of the pre -flight Monte Carlo predictions. 

To better understand these discrepancies, a systematic model reconciliation study was undertaken to determine 
the value of the aerodynamic uncertainty terms that best aligned the simulation predictions with the observed flight 
data. These uncertainty values were adjusted iteratively until the angle of attack, angle of sideslip, and body roll 
angle from the simulation results matched the flight response. Additionally, an iterative Kalman filter was developed 
to estimate the uncertainty terms using aerodynamic forces and moments derived from flight accelerations. The 
main advantage of the Kalman filter is that it had the ability to estimate the parameters simultaneously. These two 
methods were used together to understand why the profiles differed and how the mo dels could have been adjusted so 
that the Monte Carlo predictions better captured the flight behavior. 

Figure 4 illustrates that the aerodynamic un- 
certainty terms could be adjusted so that the simu- 
lation results matched the flight data reasonably 
well. These results confirmed that the slower- 
than -predicted angle -of-attack lesponse was due 
to inaccuracies in the prediction of the pitch trim 
location. While the aerodynamic uncertainty 
model that was employed in the simulation had 
adequately captured uncertainties in the magni- 
tude of the aerodynamic coefficients (i.e., slope of 
the pitching moment curve), the estimated varia- 
tion in the pitch trim location was too low be- 
cause the Cm 0 bias term had been under pre- 
dicted. Even so, the integrator gain was sufficient 
so that the controller was able to completely le- 
move the angle -of-attack error prior to the start of 
the engine test sequence. 

Despite this result, the Cmo uncertainty levels 
were not changed for the Mach 10 simulation 
analysis because CFD and wind tunnel results did 
not conclusively show that a similar behavior 
would occur at the higher Mach number. Instead, 
to allow for the possibility of similar underprediction of Cm 0 bias at Mach 10, the controller was modified by in- 
creasing the integrator gain for angle -of-attack. A stress test was conducted using the updated simulation with a con- 
servative Cm 0 uncertainty range to show that the system was robust to the differences that were observed in Mach 7 
flight. The results of this stress test are discussed in a later section. 

Another result of the aerodynamic model reconciliation study came when using the Kalman filter to estimate the 
free flight aerodynamic uncertainty parameters. Recall that in the Mach 7 simulation aerodynamic interference ef- 
fects were assumed to be completely dissipated at a horizontal separation distance of 69 inches. When the Kalman 
filter was used to estimate the free flight uncertainty model parameters to best match the aerodynamic force and 
moment flight data just outside this region, at an Xsep value of -75 inches, the filter was unable to converge. This 
process was repeated for a range of Xsep values, and the results for the axial force are shown in Fig. 5. The plot 
shows that the filter began to converge with the flight data around X S ep = -110 inches, suggesting that aerodynamic 
interference effects actually persisted until then. Consequently, the region of transition from aerodynamic interfer- 
ence to free flight was increased in the Mach 10 simulation; that is, the point at which the free flight aerodynamic 
database was used alone was extended to an X S ep value of -110 inches. In addition, a 3- a uncertainty value of +10 
inches was included for this parameter in Monte Carlo analyses. 



Figure 4. Comparison of Mach 7 RV angle- of- attack time 
history. Flight data vs. pre-flight Monte Carlo. 
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10Q 200 300 400 600 600 700 800 900 
Horizontal separation distance, in. 

Figure 5. Axial aerodynamic force comparison. Flight data 
vs. reconciled model prediction. 


The aerodynamic interference database 
(used exclusively in the first 44 inches of hori- 
zontal separation distance) was predominantly 
constructed using wind tunnel measurements 
taken at the Arnold Engineering Development 
Center at a Mach 6 test condition. This data 
was used in the Mach 7 pre -flight analysis 
without a Mach 6-to-Mach 7 adjustment. For 
the Mach 10 flight, the longitudinal coefficients 
Ca, C n and C m were adjusted to account for the 
difference in Mach by adding a constant bias 
that was based on Mach 10 CFD solutions. The 
increments that were applied were ACa = 

+0.0023, AC n = +0.0062 and AC m = -0.00075. 

Figure 6 shows the nose-down effect of the 
adjustment on the pitching moment coefficient 
in the first 44 inches of Xj EP for a nominal 
separation trajectory. In the simulation, these 
biases were linearly ramped to zero between 
Xsep = -44 and -110 in. since the free flight 
aerodynamic database was Mach dependent, 
and experimental (wind tunnel) data existed for 
Mach 6 and Mach 10 conditions. Also shown in 
the figure are the uncertainty bounds, which 
were expanded at X SEP = -44 to better capture 
the uncertainty in the interference region aero- 
dynamics. Similar increases were made in the 
uncertainty bounds for Ca and Cn- 

III. Analysis Results & Discussion 

A. Redesign of Separation Control Logic 

Control of the research vehicle during sepa- 
ration was primarily concerned with managing 
the separation upset to avoid HXLV re -contact 
and large attitude excursions that would make it 
difficult to achieve the target attitude for the 
engine test. Because of the short 2.5 sec dura- 
tion of the separation event, the effectiveness of 
applying feedback control to improve separa- 
tion performance is limited. Furthermore, the 
X-43A control system was designed primarily 

to optimize vehicle performance during the engine test; thus, it was generally left only as a contingency option to 
change the gains and structure of the controller to improve separation performance. Consequently, for the Mach 7 
flight, the separation control logic was largely a sequencer, defining when and how the various feedback loops of the 
control system were closed. In addition, there was also an elevator feed forward command. Considering the success 
of the Mach 7 separation as well as the constrained project timelines, this same structure of the separation control 
logic was retained for the Mach 10 flight. 

The Mach 10 separation control logic structure consisted of two basic elements, namely loop closure times, and 
an elevator feed-forward command. Loop closure times specify when each of the feedback loops are closed includ- 
ing roll rate, yaw rate, pitch rate, angle of attack, and roll angle feedbacks. The yaw rate loop closure is unique in 
that it can be faded in over time. This fade is specified by two times, namely fade-start and fade-end. For the Mach 7 
flight, the yaw rate fade was initiated 400 msec after separation and ended at 1000 msec. The elevator feed-forward 
command is modeled as a linear function of separation angle of attack (d e ff = bias + slope*a sep ). For the Mach 7 
flight, the bias was set to 6 degrees and the slope was 0. Although the Mach 7 structure was adopted, the flexibility 



^sepi irr 

Figure 6. Summary of changes to Mach 
interference model for Mach 10 analysis. 


7 aerodynamic 
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did exist to adjust logic values. These adjustments were important in achieving predicted separation performance on 
par with the Mach 7 pre- flight predictions. 

Re-design of the separation control logic was performed by parameterizing the structure described above and 
conducting a series of parametric studies and Monte Carlo analyses. The parameterization of the separation control 
logic resulted in the 8 design parameters: pitch rate loop (q-loop) closure time, roll rate loop (p-loop) closure time, 
yaw rate fade start time (r-fade start), yaw rate fade end time (r-fade end), angle -of -attack loop (a-loop) closure 
time, roll angle loop ((|)-loop) closure time, elevator feed-forward bias and elevator feed-forward slope. 

A systematic approach was taken to determine the best combination of the above 8 parameters. Initially a 
benchmark Monte Carlo study was performed to establish a starting point and assess the significance of the problem 
Next, parametric studies were conducted to identify promising regions of the trade space. At most, two parameters 
were varied at a time. Contour plots were used to identify regions of interest. Since these plots are based on a single 
test case (nominal or off-nominal), Monte Carlo analysis was used to assess statistical performance. Together, with 
the parametric study results it was possible to bound the potential values of each parameter. With ranges defined, 
additional time was afforded to explore candidate combinations. Again, more extensive parametric and Monte Carlo 
studies were used to ultimately select the nominal separation control logic settings. 

1 . Mach 1 0 Benchmark Monte Carlo 

The benchmark Monte Carlo study consisted of 400 perturbed Mach 10 separation cases using the Mach 7 sepa- 
ration control logic parameters. These parameters included a 6 degree elevator feed-forward bias, zero slope, q-loop 
and p-loop closures at 100 msec, a-loop and (|)-loop closures at 500 msec, a yaw rate fade-start at 400 msec, and a 
yaw rate fade-end at 1000 msec. Results of the study are shown in Table 2 and include a comparison to the final 
Mach 7 pre -flight Monte Carlo study. 


Table 2. Mach 10 Benchmark Monte Carlo Results 



Mach 10 Benchmark 

Mach 7 Monte Carlo 

Failures, % 

6.3 

2.4 

Re-contacts, % 

0.3 

0.9 

a-P Target Box, % 

89.8 

98.4 

Wing Rate Limiting, % 

27.7 

19.5 

Rudder Rate Limiting, % 

8.3 

1.5 

Wing Position Limiting, % 

< 1 

< 1 

Rudder Position Limiting, % 

6.9 

< 1 


The benchmark results show an overall degradation in separation performance at the conditions for the Mach 10 
flight. With the Mach 7 control logic settings, the probability of being within the primary a-(3 target box 
(0° = a = 1°, -0.5° = P = 0.5°) decreased, rudder rate and position limiting increased, and there was a higher overall 
failure rate. A case was considered failed if there was a re -contact, a or (3 excursions exceeded +10 deg, or the RV 
attitude at the end of separation was outside minimum acceptable ranges (0° = a = 3°, -3° = P = 3?). The only im- 
provement with respect to the Mach 7 baseline was a decrease in the probability of a re -contact which occurred be- 
cause of the favorable nose-down AC m increment derived from CFD analysis and applied in the interference region 
to account for flight at Mach 10. 

Some of the reduction in performance can be attributed to the effect of the increased Mach number on the vehi- 
cle aerodynamics. For instance, at Mach 10 the RV had less free flight stability margin and control effectiveness 
than at Mach 7. There were also performance losses due to model changes that were made in response to Mach 7 
flight data reduction, such as the increased uncertainty bands around the aerodynamic interference coefficients and 
the extension of the aerodynamic interference region, which increased the time that the vehicle was affected by 
aerodynamic interference by roughly 150 ms. Finally, changes in the RV and HXLV mass properties also had an 
effect on separation performance. 

The objective of the separation control logic redesign was to recover as much of the lost separation performance 
as possible. As will be shown, much of this performance can be recovered by judicious selection of the separation 
control logic parameters. The failure rate can be reduced by adjusting the elevator bias. Additional reductions were 
achieved through adjustment of other separation control logic parameters. To a lesser extent, the rudder rate and 
position limiting could be reduced, but not to Flight 2 levels. Some of the modeling differences between the Mach 7 
and 10 simulations contribute to this discrepancy. 
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2. Parametric Studies 

Parametric studies were performed for 
two variables at a time. For the most part, a 
nominal case was used for parametrics 
involving longitudinal parameters (elevator 
feed-forward, a- loop, and q-loop). An off- 
nominal case was lequired to exercise the 
lateral-directional dynamics of the system 
and thus enable evaluation of the p-loop, <\ > 
loop, and yaw rate fade times. 

A key study involved variation of the 
elevator feed-forward bias along with the 
a- loop closure. The other 6 control logic 
parameters were held at the benchmark 
values. Figure 7 is a typical contour plot 
produced in the parametric studies. For this 
case the color contours represent the rate 
change of angle of attack of the RV at 2.5 
seconds and the black lines represent con- 
tours of constant angle of attack at 2.5 sec- 
onds. The ideal combination is 0 rate and a 
value of 1.0 degree. An a-loop in the 
neighborhood of 500 msec provided low 
rates for entire elevator bias range. In order 
to refine the bounds on the bias, additional 
metrics were examined. 

Two of the fundamental considerations 
in selecting the elevator bias were its im- 
pact on RV/HXLV proximity and wing 
rate limiting. Proximity is defined as the 
closest distance between the RV and 
HXLV at any point in the separation trajec- 
tory and is a measure of the re -contact po- 
tential. Higher values of elevator bias force 
larger pitch down moments on the RV 
which is good from a proximity standpoint. 
However, these larger moments cause ad- 
ditional stress on the vehicle control sys- 
tem in terms of its recovery. This is k- 
flected in higher wing rates. On the other 
hand, if the elevator bias is too low, large 
wing rates are required to arrest the result- 
ing pitch up response. Ultimately the 
choice of the bias must balance these con- 
siderations. The effect of the bias value on 



300 400 500 600 700 

u-loop closure, msec 

Figure 7. Parametric study results to determine acceptable 
elevator bias and a-loop closure ranges. 
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400 500 600 700 800 900 1000 

a-loop closure, msec 

Figure 8. Parametric study results to determine effect of elevator 
bias and a-loop closure on rudder rate limiting. 


the wing rate is shown in Fig. 8 where the black lines on the contour plot define the range of elevator bias (roughly 
2.5 to 5.5 deg) that limits peak wing rates to +50 deg/s. Figure 9 shows the relationship between elevator bias and 
proximity. A spread in proximity of roughly 1 inch can be seen across the range of biases that were exa mined. 

Another key parametric study was conducted to determine the yaw rate fade parameters. The tendency of the RV 
to rely on considerable rudder deflection in response to yawing motion has long been recognized. 6 The large aero- 
dynamic uncertainty associated with separation interference makes rudder activity particularly high during the sepa- 
ration event. In order to avoid the potential complications of rudder rate and position limiting the separation control 
logic allows the yaw rate feedback to be faded in over time. Effectively, design of the fade reduces to a trade be- 
tween rudder activity and sideslip/roll angle excursions. 
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Figure 9. Effect of elevator bias setting on proximity at 
axial horizontal distance of 37 inches. 


For this study, an off-nominal case 
was chosen that had large rudder 
deflections and a peak (3 of over 5 deg. 

This same case was run with a range 
of fade start times and duration. Figure 
10 illustrates the results where color 
contours reflect maximum rudder 
position and the black traces are 
associated with contours of peak 
sideslip angle. Two regions of interest 
that have relatively low rudder activity 
are shaded. The first is associated with 
early start times and fade durations in 
the range of 300 to 800 msec. The 
second region is associated with later 
start times and mo derate to high fade 
durations. The peak sideslip angle, and 
similarly roll angle (not pictured), 
varied predominantly with the fade 
start time, and was farily insensitive to 
the fade duration, Thus, lower fade 
start times led to lower values of (3 max . 

This explains why the fade start times 
associated with the second egion is 
limited to 500 to 600 msec. These two 
regions became the focus of additional 
Monte Carlo studies described below. 

3. Monte Carlo Studies 

Monte Carlo studies were used to 
compare the candidate settings related 
to the yaw rate fade. The two regions 
of interest previously identified in the 
yaw rate fade parametric study were 
targeted. Over 40 separate Monte 
Carlo analyses were run at a range of 
fade-start and duration settings within 
these regions. In addition, variation in 
the (|)-loop (100 and 500 msec) and 
elevator feed-forward bias (1 and 4 
deg.) was included to capture potential 
coupling between logic parameters. 

The results showed that the two le- 
gions represented a choice between 
emphasizing lower failure rate vs. 
lower control surface activity. The 

results of a representative case from each region are shown in Table 3. By shifting the fade later in the separation 
timeline, the rudder rate and deflection limiting could be significantly decreased, but only with an accompanying 
degradation in the mission success statistics. 

Monte Carlo studies were also used to quickly understand whether there was any advantage to delaying the clo- 
sure of pitch and roll rate loops beyond the benchmark value of 100 msec. For both of these parameters, loop closure 
times less than 100 msec were not considered because the RV does not leave the separation pistons until roughly 
that time. Two separate Monte Carlo studies were performed for both loops with closure times of 200 msec. In both 
cases separation performance degraded as loop closure times increased beyond 100 msec. Further Monte Carlo 
analysis was not performed since delaying rate feedback control any longer would only tend to further degrade per- 
formance. 



400 500 600 700 
Fade staff msec 

Figure 10. Parametric study results to identify promising regions of 
yaw fade parameter trade space. 
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Table 3. Statistical Performance for Two Candidate Yaw Fade Strategies 


Metric 

Elev. bias = 3.5° 
Fade start = 100 ms 
Fade end = 700 ms 
(400 cases) 

Elev. bias = 3.5° 
Fade start = 600 ms 
Fade end = 1500 ms 
(400 cases) 

Failure Rate 

3.5% 

5.5% 

Within Primary a-(3 Target Box 

95.5% 

89.5% 

iPmaxI ^ 4 deg 

5.2% 

8.2% 

Wing Rate >100 deg/s 

23.1% 

21.7% 

Rudder Rate >100 deg/s 

13.5% 

1.3% 

Reach Rudder Position Limits 

4.9% 

1.6% 


4. Mission Data Load 

Based on the results of these studies, a range of discrete values were defined for each of these parameters and in- 
cluded in the Mission Data Load (MDL). With the MDL the capability existed to upload selected flight control, 
navigation and propulsion parameters to the RV at mission time instead of early in the program when the flight 
software was released. The advantage of this arrangement is that parameters included in the MDL could be finalized 
much later in the project timeline; however, more MDL parameters meant that additional verification tests had to be 
run. For the Mach 10 flight, the separation logic parameters were added to the MDL, effectively removing further 
separation logic re-design from the critical path. However, the range of potential values had to be chosen carefully to 
limit the extent of verification testing needed to qualify the Mission Data Load of the Flight Management Unit. The 
ranges and nominal values are specified in Table 4. The nominal slope and a-loop closure times were selected from 
parametric study results. Monte Carlo comparisons were used to select the elevator bias, (|>loop and yaw fade times. 
The statistical results of each Monte Carlo were combined in a cost function to determine which set of parameters 
performed best. The cost function was weighted to account for different objectives, with the final selection made 
using a cost function that weighted failure rate and target box accuracy five times more than surface rate and deflec- 
tion limiting. 


Table 4. Mission Data Load Ranges of Separation Control Parameters. Nominal Values in Red. 
Variable Values 


Bias, deg 

1.0 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 

5.5 

Slope 

-2 

-1 

0 

1 







q-loop closure , msec 

100 










a-loop closure, msec 

100 

300 

500 

700 







p-loop closure, msec 

100 










(|)-loop closure, msec 

100 

300 

500 

700 







Start yaw- loop fade, msec 

100 

200 

400 

500 

600 






End yaw- loop fade, msec 

500 

600 

700 

800 

900 

1000 

1300 

1400 

1500 

1600 


B. Mach 10 Preflight Monte Carlo Analysis 

Table 5 compares the key results from the final pre- flight Mach 10 Monte Carlo analysis with the pre -flight 
Mach 7 values. Both sets of statistical results were derived from 1200 off-nominal cases. The results show an im- 
provement over the Mach 10 benchmark case and demonstrate that much of the lost separation performance was 
recovered by adjusting the separation control logic settings. As expected, there were more failures and lower target 
box accuracy for the Mach 10 results when compared to the Mach 7 results. The number of failed cases increased 
primarily because there were more cases with a and (3 excursions exceeding 10 deg. The biggest cause of the in- 
crease in large-excursion cases was the longer interference-to-free transition region since it required more extrapola- 
tion of the 8-dimensional interference aerodynamic database. Also, there were actually less re -contacts for the Mach 
10 case, due primarily to the effect of the Mach 7-to-Mach 10 increments applied in the aerodynamic interference 
region which tended to increase the upward normal force and induce a more nose-down pitching moment on the RV. 
While the frequency of wing rate limiting increased only moderately in the Mach 10 results, there was significantly 
more rudder rate and position limiting than in the Mach 7 results. Some of this additional rudder activity was also 
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Table 5. Final Mach 10 Pre-Flight Monte Carlo Results 



Mach 10 Baseline 

Mach 7 Monte Carlo 

Number of Cases 

1200 

1200 

Failures, % 

4.4 

2.2 

Re-contacts, % 

0.2 

0.9 

Primary Target Box, % 

95.0 

98.3 

Minimum a (average), deg 

-1.64 

-1.40 

Peak P (average +/-), deg 

1.12/-0.80 

0.94/-0.77 

Wing Rate Limiting, % 

24.5 

19.5 

Rudder Rate Limiting, % 

8.6 

1.5 

Wing Position Limiting, % 

< 1 

< 1 

Rudder Position Limiting, % 

4.4 

< 1 


caused by the increased extrapolation of interfer- 
ence aerodynamic coefficients brought on by the 
longer transition region. 

C. Cmo bias Stress Test 

A stress test was conducted with the Mach 10 
simulation to ensure that the system was robust to 
the type of Cmo bias error that was observed in the 
Mach 7 flight. First, two individual cases were run 
with the free flight pitching moment increment set 
to +0.002; which was larger than the value that was 
deduced from the Mach 7 model reconciliation 
study. The resulting angle -of-attack profiles are 
shown in Fig. 1 1 and indicate that the RV is capable 
of recovering from either extreme in pitching mo- 
ment bias. In addition, a Monte Carlo analysis was 
performed in which the free flight Cmo increment 
was varied uniformly between -0.002 and +0.002, 
and all other variables were varied using baseline 
uncertainty ranges. The results were similar to the 
results of the final pre -flight Monte Carlo, with only 
2 additional failures in the CmO set. The resulting 
Monte Carlo uncertainty envelope for angle -of- 
attack is also shown in Fig. 1 1 . 

IV. Mach 10 Flight Performance 

On Nov. 16, 2004 the final X43A flight oc- 
curred with a fully successful stage separation. 11 " 12 
The HXLV accurately delivered the RV to the sepa- 
ration conditions well within one standard deviation 
of the target values. The RV exhibited nearly nomi- 
nal performance during separation and based on 
analysis of flight data there was no indication of a 
re -contact. Figure 12 shows the angle -of-attack pro- 
file reconstructed from flight telemetry as well as 
the predicted profile from the simulation (with 
flight-derived initial conditions). As in the Mach 7 
flight, the RV took longer than predicted to reach 
the commanded angle -of-attack, again, most likely 



Figure 11. Results of Mach 10 Cmo bias stress test. 



Figure 12. Comparison of Mach 10 RV angle-of- attack 
time history. Flight data vs. pre- flight Monte Carlo. 
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due to inaccuracies in the prediction of the pitch 
trim location. Even so, the RV still achieved the 
required angle -of-attack prior to the engine test, and 
did so quicker than in the Mach 7 flight because of 
the controller modification that increased the angle - 
of-attack integrator gain. Finally, note that the entire 
angle -of-attack response was captured within the 
pre -flight Monte Carlo uncertainty bound. 

Figure 13 shows the elevator deflection during 
separation for the Mach 10 flight. Again, the flight 
data was well captured within the pre -flight Monte 
Carlo envelope. The small bias after the transient is 
due to the error in CmO, resulting from the incorrect 
trim location in the simulation aerodynamic data- 
base. There was very little lateral-directional e- 
sponse during the Mach 10 separation, indicating 
that both pistons fired nearly simultaneously. The 
sideslip angle time history and pre -flight Monte 
Carlo bounds are shown in Fig. 14. 



Figure 13. Comparison of Mach 10 RV elevator deflection 
time history. Flight data vs. pre- flight Monte Carlo. 


Nominal Prediclion_ 

Flight Data 

Pre-IUghl 
Monle Carlo 



V. Conclusion 

The success of stage separation during the first 
powered flight of the X-43A in March 2004 pro- 
vided confidence in the separation analysis ap- 
proach and simulation modeling. For the next flight 
conducted at Mach 10 conditions, a similar ap- 
proach was undertaken and the same simulation 
structure was adopted. Considerable effort was 
made to identify Mach 10 model updates through 
Mach 7 flight data reduction, model reconciliation 
studies, and CFD analysis. Once these updates were 
made, benchmark analysis of the Mach 10 mission 
showed degraded performance relative to Mach 7 
pre -flight predictions. Much of the performance was 
recovered by adjusting the feed-forward elevator 
command and feedback loop closure times. This 
redesign of the separation control logic was 
achieved through extensive parametric studies and 
Monte Carlo analyses conducted using POST. In 

addition, an independent, high-fidelity ADAMS simulation was also maintained throughout the analysis cycle to 
understand the effects of high-order models and to verify the POST simulation. Numerous Monte Carlo trade studies 
with conservative uncertainty models were performed with these tools to understand the impact of modeling as- 
sumptions and to perform stress testing to ensure that the system was robust to anomalies observed in the Mach 7 
flight. 

On November 16, 2004 the X-43A Mach 10 flight occurred with another successful stage separation. The Re- 
search Vehicle exhibited near-nominal performance throughout separation and the vehicle performed within the 
uncertainty bounds predicted in the pre -flight analysis and went on to break the world record for airbreathing pow- 
ered flight, set previously by the Mach 7 X-43A flight. 


1.0 1.5 

Time, sec 

Figure 14. Comparison of Mach 10 RV sideslip angle time 
history. Flight data vs. pre-flight Monte Carlo. 
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